TITLE OF THE INVENTION 

Radio Equipment Capable of Real Time Change of Antenna 
Directivity and Doppler Frequency Estimating Circuit Used for the Radio 
Equipment 

5 BACKGROUND OF THE INVENTION 
Field of the Invention 

The present invention relates to a radio equipment that can change 
antenna directivity on real time basis as well as to a Doppler frequency 
estimating circuit used therefor. More specifically, the present invention 
10 relates to a radio equipment used in an adaptive array radio base station, 
and to a Doppler frequency estimating circuit used therefor. 
Description of the Background Art 

Recently, various methods of transmission channel allocation have 
been proposed to realize effective use of frequency, in a mobile 
15 communication system, of which some have been practically implemented. 

Fig. 15 shows an arrangement of channels in various communication 
systems including frequency division multiple access (FDMA), time division 
multiple access (TDMA) and path division multiple access (PDMA). 

Referring to Fig. 15, FDMA, TDMA and PDMA will be briefly 
20 described. Fig. 15(A) represents FDMA in which analog signals of users 1 
to 4 are subjected to frequency division and transmitted over radio waves of 
different frequencies fl to f4, and the signals of respective users 1 to 4 are 
separated by frequency filters. 

In TDMA shown in Fig. 15(B), digitized signals of respective users 
25 are transmitted over the radio waves having different frequencies fl to f4 
and time-divided time slot by time slot (time slot: a prescribed time period), 
and the signals of respective users are separated by the frequency filters and 
time-synchronization between a base station and mobile terminals of 
respective users. 

30 Recently, PDMA method has been proposed to improve efficiency of 

use of radio frequency, as portable telephones have come to be widely used. 
In the PDMA method, one time slot of one frequency is spatially divided to 
enable transmission of data of a plurality of users, as shown in Fig. 15(C). 
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In the PDMA, signals of respective users are separated by the frequency 
filters, the time synchronization between the base station and the mobile 
terminals of respective users, and a mutual interference eliminating 
apparatus such as an adaptive array. 
5 The operation principle of such an adaptive array radio base station 

is described in the following literature, for example: 

B. Widrow, et al.: "Adaptive Antenna Systems", Proc. IEEE, vol.55, 
No.l2,pp.2143-2159 (Dec. 1967). 

S. P. Applebaum: "Adaptive Arrays", IEEE Trans. Antennas & 
10 Propag., vol.AP-24, No.5, pp.585-598 (Sept. 1976). 

O. L. Frost, III: " Adaptive Least Squares Optimization Subject to 
Linear Equality Constraints", SEL- 70-055, Technical Report, No.6796-2, 
Information System Lab., Stanford Univ. (Aug. 1970). 

B. Widrow and S. D. Stearns: "Adaptive Signal Processing", 
15 Prentice-Hall, Englewood Cliffs (1985). 

R. A. Monzingo and T. W. Miller: "Introduction to Adaptive Arrays", 
John Wiley & Sons, New York (1980). 

J. E. Hudson: "Adaptive Array Principles", Peter Peregrinus Ltd., 
London (1981). 

20 R. T. Compton, Jr.: "Adaptive Antennas — Concepts and 

Performance", Prentice-Hall, Englewood Cliffs (1988). 

E. Nicolau and D. Zaharia: "Adaptive Arrays", Elsevier, 
Amsterdam (1989). 

Fig. 16 is a model diagram conceptually showing the operation 
25 principle of such an adaptive array radio base station. Referring to Fig. 16, 
an adaptive array radio base station 1 includes an array antenna 2 formed 
by n antennas #1, #2, #3, -, #n, and a first diagonal line area 3 shows a 
range in which radio waves from the array antenna 2 can be received. A 
second diagonal line area 7 shows a range in which radio waves from 
30 adjacent another radio base station 6 can be received. 

In the area 3, the adaptive array radio base station 1 
transmits/receives a radio signal to/from a portable telephone 4 forming a 
terminal of a user A (arrow 5). In the area 7, the radio base station 6 
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transmits/receives a radio signal to/from a portable telephone 8 forming a 
terminal of another user B (arrow 9). 

When the radio signal for the portable telephone 4 of the user A 
happens to be equal in frequency to the radio signal for the portable 
5 telephone 8 of the user B, it follows that the radio signal from the portable 
telephone 8 of the user B serves as an unnecessary interference signal in the 
area 3 depending on the position of the user B, to disadvantageous^ mix 
into the radio signal transmitted between the portable telephone 4 of the 
user A and the adaptive array radio base station 1. 

10 In this case, it follows that the adaptive array radio base station 1 

receiving the mixed radio signals from both users A and B in the 
aforementioned manner outputs the signals from the users A and B in a 
mixed state unless some necessary processing is performed, to 
disadvantageously hinder communication with the regular user A. 

15 [Configuration and Operation of Conventional Adaptive Array 

Antenna] 

In order to eliminate the signal from the user B from the output 
signal, the adaptive array radio base station 1 performs the following 
processing. Fig. 17 is a schematic block diagram showing the structure of 
20 the adaptive array radio base station 1. 

Assuming that A(t) represents the signal from the user A and B(t) 
represents the signal from the user B, a signal xl(t) received in the first 
antenna #1 forming the array antenna 2 shown in Fig. 16 is expressed as 
follows: 

25 xl(t) = al x A(t) + bl x B(t) 

where al and bl represent factors changing in real time, as described later. 
A signal x2(t) received in the second antenna #2 is expressed as 

follows: 

x2(t) = a2 x A(t) + b2 x B(t) 
30 where a2 and b2 also represent factors changing in real time. 

A signal x3(t) received in the third antenna #3 is expressed as 

follows: 

x3(t) = a3 x A(t) + b3 x B(t) 
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where a3 and b3 also represent factors changing in real time. 

Similarly, a signal xn(t) received in the n-th antenna #n is expressed 
as follows: 

xn(t) = an x A(t) + bn x B(t) 
5 where an and bn also represent factors changing in real time. 

The above factors al, a2, a3, -, an show that the antennas #1, #2, #3, 
-, #n forming the array antenna 2 are different in receiving strength from 
each other with respect to the radio signal from the user A since the relative 
positions of the antennas #1, #2, #3, #n are different from each other (the 
10 antennas #1, #2, #3, #n are arranged at intervals about five times the 
wavelength of the radio signal, i.e., about 1 m, from eath other). 

The above factors bl, b2, b3, •-, bn also show that the antennas #1, 
#2, #3, #n are different in receiving strength from each other with respect 
to the radio signal from the user B. The users A and B are moving and 
15 hence these factors al, a2, a3, an and bl, b2, b3, bn change in real 
time. 

The signals xl(t), x2(t), x3(t), xn(t) received in the respective 
antennas #1, #2, #3, #n are input to a receiving unit 1R forming the 
adaptive array radio base station 1 through corresponding switches 10-1, 

20 10-2, 10-3, 10-n respectively so that the received signals are supplied to a 
weight vector control unit 11 and to one inputs of corresponding multipliers 
12-1, 12-2, 12-3, 12-n respectively. 

Weights wl, w2, w3, wn for the signals xl(t), x2(t), x3(t), -, xn(t) 
received in the antennas #1, #2, #3, #n are applied from the weight vector 

25 control unit 11 to other inputs of these multipliers 12-1, 12-2, 12-3, •-, 12-n 
respectively. The weight vector control unit 1 1 calculates these weights wl, 
w2, w3, wn in real time, as described later. 

Therefore, the signal xl(t) received in the antenna #1 is converted to 
wl x (alA(t) + blB(t)) through the multiplier 12-1, the signal x2(t) received 

30 in the antenna #2 is converted to w2 x (a2A(t) + b2B(t)) through the 

multiplier 12-2, the signal x3(t) received in the antenna #3 is converted to 
w3 x (a3A(t) + b3B(t)) through the multiplier 12-3, and the signal xn(t) 
received in the antenna #n is converted to wn x (anA(t) + bnB(t)) through the 
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multiplier 12-n. 

An adder 13 adds the outputs of these multipliers 12-1, 12-2, 12-3, -, 
12-n, and outputs the following signal: 

wl(alA(t) + blB(t)) + w2(a2A(t) + b2B(t)) + w3(a3A(t) + b3B(t)) + ... + 
wn(anA(t) + bnB(t)) 

This expression is classified into terms related to the signals A(t) and 
B(t) respectively as follows: 

(wlal + w2a2 + w 3a3 + - + wnan)A(t) + (wlbl + w2b2 + w3b3 + - + 
wnbn)B(t) 

As described later, the adaptive array radio base station 1 identifies 
the users A and B and calculates the aforementioned weights wl, w2, w3, 
wn to be capable of extracting only the signal from the desired user. 
Referring to Fig. 17, for example, the weight vector control unit 11 regards 
the factors al, a2, a3, an andbl, b2, b3, • •-, bn as constants and calculates 
the weights wl, w2, w3, wn so that the factors of the signals A(t) and B(t) 
are 1 and 0 as a whole respectively, in order to extract only the signal A(t) 
from the intended user A for communication. 

In other words, the weight vector control unit 11 solves the following 
simultaneous linear equations, thereby calculating the weights wl, w2, w3, 
wn on real time basis so that the factors of the signals A(t) and B(t) are 1 
and 0 respectively: 

wlal + w2a2 + w3a3 + •- + wnan = 1 

wlbl + w2b2 + w3b3 + •• + wnbn = 0 

The method of solving the above simultaneous linear equations, not 
described in this specification, is known as described in the aforementioned 
literature and already put into practice in an actual adaptive array radio 
base station. 

When setting the weights wl, w2, w3, wn in the aforementioned 
manner, the adder 13 outputs the following signal: 
output signal = 1 x A(t) + 0 x B(t) = A(t) 
[User Identification, Training Signal] 

The aforementioned users A and B are identified as follows: 
Fig. 18 is a schematic diagram showing the frame structure of a 



radio signal for a portable telephone set. The radio signal for the portable 
telephone set is roughly formed by a preamble consisting of a signal 
sequence known to the radio base station and data (sound etc.) consisting of 
a signal sequence unknown to the radio base station. 
5 The signal sequence of the preamble includes a signal sequence of 

information for recognizing whether or not the user is a desired user for 
making communication with the radio base station. The weight vector 
control unit 11 (Fig. 17) of the adaptive array radio base station 1 compares 
a training signal corresponding to the user A fetched from a memory 14 with 

10 the received signal sequence and performs weight vector control (decision of 
weights) for extracting a signal apparently including the signal sequence 
corresponding to the user A. The adaptive array radio base station 1 
outputs the signal from the user A extracted in the aforementioned manner 
as an output signal Sitx(t). 

15 Referring again to Fig. 17, an external input signal STx(t) is input to 

a transmission unit IT forming the adaptive array radio base station 1 and 
supplied to one inputs of multipliers 15-1, 15-2, 15-3, 15-n. The weights 
wl, w2, w3, wn previously calculated by the weight vector control unit 11 
on the basis of the received signal are copied and applied to other inputs of 

20 these multipliers 15-1, 15-2, 15-3, 15-n respectively. 

The input signals Stx(0 weighted by these multipliers 15-1, 15-2, 
15-3, •-, 15-n are sent to the corresponding antennas #1, #2, #3, •», #n 
through corresponding switches 10-1, 10-2, 10-3, 10-n respectively, and 
transmitted into the area 3 shown in Fig. 16. 

25 The signal transmitted through the same array antenna 2 as that in 

receiving is weighted for the target user A similarly to the received signal, 
and hence the portable telephone set 4 of the user A receives the transmitted 
radio signal as if the signal has directivity to the user A. Fig. 19 images 
such transfer of a radio signal between the user A and the adaptive array 

30 radio base station 1. Imaged is such a state that the adaptive array radio 
base station 1 transmits the radio signal with directivity toward the target 
portable telephone set 4 of the user A as shown in a virtual area 3a in Fig. 19 
in contrast with the area 3 of Fig. 16 showing the range actually receiving 



radio waves. 

As described above, in the PDMA method, a technique is necessary 
to remove co-channel interference. In this point, an adaptive array that 
places nulls on the interfering waves adaptively is an effective means, as it 
5 can effectively suppress the interfering wave even when the level of the 
interfering wave is higher than the level of the desired wave. 

When an adaptive array is used for a base station, it becomes 
possible not only to remove interference at the time of reception but also to 
reduce unnecessary radiation at the time of transmission. 

10 At this time, an array pattern at the time of transmission may be an 

array pattern for reception, or the array pattern may be newly generated 
based on a result of incoming direction estimation, for example. The latter 
method is applicable no matter whether FDD (Frequency Division Duplex) 
or TDD (Time Division Duplex) is used. It requires, however, a complicated 

15 process. When the former approach is to be used while FDD is utilized, 

modification of the array arrangement or weight becomes necessary, as the 
array patterns for transmission and reception are different. Therefore, 
generally, application is on the premise that TDD is utilized, and in an 
environment where external slots are continuous, satisfactory characteristic 

20 has been ensured. 

As described above, in the TDD/PDMA method using an adaptive 
array in the base station, when an array pattern (weight vector pattern) 
obtained for the up link is used for the down link, transmission directivity 
may possibly be degraded in the down link because of time difference 

25 between the up and down links, assuming a dynamic Rayleigh propagation 
degree with angular spread. 

More specifically, there is a time interval from transmission of the 
radiowave from a user terminal to the base station through the up link until 
radiowave is emitted from the base station to the user terminal through the 

30 down link. Therefore, if the speed of movement of the user terminal is not 
negligible, transmission directivity degrades because of the difference 
between the direction of radiowave emission from the base station and the 
actual direction of the user terminal. 
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As a method of estimating weight for the down link considering such 
a variation in the propagation path, a method of estimating a transmission 
response vector of the down link by extrapolation utilizing a reception 
response vector obtained in the up link has been proposed. 
5 When there is an estimation error in the reception response vector 

estimated for the up link because of a noise in the reception signal or 
sampling error, there would be an error in the result of extrapolation, 
making it difficult to correctly estimate the transmission response vector for 
the down link, and hence, it becomes impossible to realize satisfactory 

10 control of transmission directivity. 

Now, propagation environment of the propagation path is 
represented by variation of reception coefficient of the propagation path, 
that is, degree of fading. The degree of fading is represented by a so-called 
Doppler frequency (FD) as a physical amount. 

15 More specifically, dependent on the degree of fading on the 

propagation path, there would be an error in the result of extrapolation. 
Therefore, in order to prevent extrapolation error, it is necessary to know the 
degree of fading of the propagation path, that is, to know the Doppler 
frequency. 

20 A method has been proposed as disclosed, for example, in Japanese 

Patent Laying-Open No. 7-162360, in which the degree of fading is 
estimated by finding a correlation value of reference signals included in 
reception signals received preceding and succeeding in time with respect to 
each other. In the conventional method, however, the correlation value is 

25 calculated using reference signals that are included in the reception signals 
themselves, which means that there are much interfering components, 
making correct estimation difficult. 

Further, the timing of the reference signals is fixed. Therefore, it 
was impossible to calculate the correlation value at an arbitrarily timing, 

30 and therefore, the calculating process was not flexible. 

Meanwhile, the method of estimating Doppler frequency of the 
propagation path for each user terminal separated by the adaptive array 
processing has not yet been developed. 
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SUMMAKY OF THE INVENTION 

An object of the present invention is to provide a radio equipment 
capable of correctly estimating a transmission response vector of the down 
link even when there is an estimation error in the reception response vector 
5 estimated for the up link, and hence capable of satisfactory control of 
transmission directivity. 

Another object of the present invention is to provide a Doppler 
frequency estimating circuit for estimating the Doppler frequency for each 
user terminal, which is free from the influence of the interfering component 

10 of the reception signal and has increased flexibility in an operating process. 

A still further object of the present invention is to provide a radio 
equipment employing the Doppler frequency estimating circuit for 
estimating the Doppler frequency of each user terminal, which is free from 
the influence of the interfering component of the reception signal and has 

15 increased flexibility in an operating process. 

According to an aspect, the present invention provides a radio 
equipment transmitting/receiving signals time-divisionally to/from a 
plurality of terminals while changing antenna directivity on real time 
basis, including a plurality of antennas arranged in a discrete manner, and 

20 a transmission circuit and a reception circuit that share the plurality of 

antennas for transmitting/receiving signals. The reception circuit includes 
a reception signal separating unit separating a signal from a specific 
terminal among the plurality of terminals, based on signals from the 
plurality of antennas when the reception signal is received, and a reception 

25 propagation path estimating unit estimating a reception response vector of a 
propagation path from the specific terminal, based on the signals from the 
plurality of antennas when the reception signal is received. The 
transmission circuit includes a transmission propagation path estimating 
unit estimating a transmission response vector of a transmission path when 

30 a transmission signal is transmitted, based on the result of estimation by the 
reception propagation path estimating unit, and a transmission directivity 
control unit updating the antenna directivity when the transmission signal 
is transmitted, based on a result of estimation by the transmission 
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propagation path estimating unit. The transmission propagation path 
estimating unit has an extrapolation processing unit calculating the 
transmission response vector of a down link slot to the specific terminal, by 
an extrapolation process based on a plurality of reception response vectors of 
5 up link slots from the specific terminal estimated by the reception 
propagation path estimating unit, a memory holding a plurality of 
parameters used for the extrapolation process, determined in advance in 
accordance with the propagation environment of the propagation path, and 
a selecting unit estimating the propagation environment of the propagation 
10 path, selecting a parameter corresponding to the estimated propagation 
environment among the held plurality of parameters, and applying the 
selected parameter to the extrapolation process by the extrapolation 
processing unit. 

Therefore, even when there is an estimation error in the reception 

15 response vector estimated for the up link, the transmission response vector 
for the down link can be estimated correctly and hence satisfactory control of 
transmission directivity can be realized, as the parameter used for the 
extrapolation process is selected in accordance with the propagation 
environment of the propagation path. 

20 Preferably, the parameter is an extrapolation distance in the 

extrapolation process by the extrapolation processing unit. The memory 
holds a plurality of extrapolation distances determined in advance in 
accordance with the Doppler frequencies representing the propagation 
environment. The selecting unit estimates the Doppler frequency of the 

25 propagation path, selects the extrapolation distance corresponding to the 
estimated Doppler frequency among the held plurality of extrapolation 
distances and applies the selected extrapolation distance to the 
extrapolation process by the extrapolation processing unit. 

Therefore, when there is an estimation error in the reception 

30 response vector of the up link, there would be a larger extrapolation error as 
the extrapolation distance for the extrapolation process becomes longer, and 
in view of this phenomenon, it becomes possible to estimate the transmission 
response vector correctly by selecting the extrapolation distance in 
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accordance with the Doppler frequency representing the propagation 
environment. 

More preferably, the selecting unit selects a shorter extrapolation 
distance when the estimated Doppler frequency is lower, and selects a 
5 longer extrapolation distance when it is higher. 

Therefore, as the variation in the propagation environment is 
smaller as the Doppler frequency is lower, extrapolation exceeding the 
actual amount of variation is prevented by making shorter the extrapolation 
distance, and as the variation in propagation environment is larger as the 
10 Doppler frequency is higher, sufficient extrapolation is ensured by making 
longer the extrapolation distance. 

More preferably, the parameter is an extrapolation distance in an 
extrapolation process by the extrapolation processing unit. The memory 
holds a plurality of extrapolation distances determined in advance in 
15 accordance with a signal error between the separated signals and an 

expected desired signals representing the propagation environment. The 
selecting unit estimates signal error of the propagation path, selects the 
extrapolation distance corresponding to the estimated signal error among 
the held plurality of extrapolation distances, and applies the selected 
20 extrapolation distance to the extrapolation process by the extrapolation 
processing unit. 

When there is a large signal error, the estimation error of the 
reception response vector of the up link becomes larger and the 
extrapolation error also becomes larger, and in view of this phenomenon, an 
25 extrapolation distance in accordance with the signal error representing the 
propagation environment is selected, whereby it becomes possible to 
correctly estimate the transmission response vector. 

More preferably, the selecting unit selects a shorter extrapolation 
distance when the estimated signal error is larger, and selects a larger 
30 extrapolation distance when it is smaller. 

When the signal error is larger, the extrapolation error increases, 
and therefore, extrapolation error is suppressed by making shorter the 
extrapolation distance, and when the signal error becomes smaller, the 
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extrapolation error becomes smaller, and therefore, sufficient extrapolation 
can be done by making longer the extrapolation distance. 

More preferably, the parameter is an extrapolation distance in an 
extrapolation process by the extrapolation processing unit. The memory 
5 holds a plurality of extrapolation distances determined in advance in 
accordance with Doppler frequencies and a signal error between the 
separated signal and an expected desired signal, which represent the 
propagation environment. The selecting unit estimates the Doppler 
frequency and the signal error of the propagation path, selects an 
10 extrapolation distance corresponding to the estimated Doppler frequency 

and the signal error among the held plurality of extrapolation distances, and 
applies the selected one to the extrapolation process by the extrapolation 
processing unit. 

Therefore, in view of the fact that the Doppler frequency and the 
15 signal error increase extrapolation error, an extrapolation distance in 

accordance with the Doppler frequency and the signal error representing the 
propagation environment is selected, whereby correct estimation of 
transmission response vector becomes possible. 

More preferably, the selecting unit temporarily selects an 
20 extrapolation distance corresponding to the estimated Doppler frequency, 
and corrects the temporarily selected extrapolation distance in accordance 
with the estimated signal error. 

Therefore, basic extrapolation distance is selected in accordance 
with the Doppler frequency that has much influence on the extrapolation 
25 error, and the extrapolation distance is corrected based on the signal error, 
whereby more accurate estimation of transmission response vector becomes 
possible. 

More preferably, the relation between the propagation environment 
and the plurality of parameters is determined individually for every radio 
30 equipment. 

Therefore, as the correspondence between the propagation 
environment and the parameters is found by measurement in advance for 
each radio equipment, more accurate estimation of transmission response 
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vector becomes possible. 

More preferably, the relation between the propagation environment 
and the plurality of parameters is determined commonly to a plurality of the 
radio equipments. 

5 Therefore, when there is only a small difference between the 

individual radio equipments, the correspondence between the propagation 
environment and the parameters is made common among a plurality of 
radio equipments, whereby the process steps for manufacturing the radio 
equipments can be simplified. 

10 According to another aspect, the present invention provides, in a 

radio equipment that changes antenna directivity on real time basis and 
transmits/receives signals time divisionally to/from a plurality of terminals, 
a Doppler frequency estimating circuit estimating Doppler frequency of a 
propagation path with a specific terminal, which includes: a reception signal 

15 separating unit separating a signal from the specific terminal among the 
plurality of terminals, based on signals received by a plurality of antennas 
arranged in a discrete manner; a reception propagation path estimating unit 
estimating a reception response vector of a propagation path from the 
specific terminal based on signals received by the plurality of antennas; a 

20 correlation operating unit calculating a vector correlation value based on 
reception response vectors preceding and succeeding in time estimated by 
the reception propagation path estimating unit; and an estimating unit 
estimating a Doppler frequency corresponding to the vector correlation 
value calculated by the correlation operating unit, based on a 

25 correspondence between vector correlation values and Doppler frequencies 
determined in advance experimentally. 

Therefore, not the reception signal itself but a correlation value of 
reception response vectors is calculated, and therefore, Doppler frequency of 
the propagation path for each separated specific terminal can be correctly 

30 estimated, free from the influence of the interfering component. 

Preferably, the correlation operating unit includes a calculating unit 
calculating an instantaneous correlation value between the reception 
response vectors preceding and succeeding in time, and outputting 
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calculated value as the vector correlation value. 

Therefore, it is possible to find correlation between reception 
response vectors at an arbitrary timing, and hence instantaneous Doppler 
frequency of the propagation path can be correctly estimated. 
5 Preferably, the correlation operating unit includes a calculating unit 

calculating an instantaneous correlation value between the reception 
response vectors preceding and succeeding in time, and an averaging unit 
weight-averaging a past correlation value and a present correlation value 
calculated by the calculating unit, with a prescribed weight coefficient, and 
10 outputting an obtained average value as said vector correlation value. 

Therefore, even when there is generated an error in the 
instantaneous Doppler frequency because of an abrupt fading, the Doppler 
frequency can be estimated correctly without the influence of such error, as 
the correlation values are averaged. 
15 More preferably, the prescribed weight coefficient is set such that a 

weight for a past correlation value is large and a weight for a present 
correlation value is small. 

As the weight of the past correlation value is made large at the time 
of averaging, even when there is generated an error in the instantaneous 
20 Doppler frequency because of an abrupt fading, more accurate estimation of 
the Doppler frequency is possible, not influenced by the error. 

More preferably, the correlation operating unit calculates a vector 
correlation value, based on a reception response vector of a present frame 
slot and a reception response vector of an immediately preceding frame slot. 
25 Therefore, different from the reference signal in the prior art, 

mutual correlation between reception response vectors can be obtained at 
any timing, and therefore, flexibility of operation to calculate the correlation 
value is increased. 

More preferably, the correlation operating unit calculates a vector 
30 correlation value based on a reception response vector of a present frame slot 
and a reception response vector of a most recent slot free of any reception 
error among past frame slots. 

Therefore, different from the reference signal in the prior art, 
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mutual correlation between reception response vectors can be obtained at 
any timing, and therefore, flexibility of operation to calculate the correlation 
value can be increased and, particularly, the influence of reception error can 
be eliminated. 

5 More preferably, the correlation operating unit calculates a vector 

correlation value based on a reception response vector of a former half and a 
reception response vector of a latter half of one slot. 

Therefore, different from the reference signal in the prior art, 
mutual correlation of the reception response vectors can be obtained at any 

10 timing, and therefore, flexibility of operation to calculate the correlation 
value can be increased. 

According to a still further aspect, the present invention provides a 
radio equipment that changes antenna directivity on a real time basis and 
transmits/receives signals time divisionally to/from a plurality of terminals, 

15 including a plurality of antennas arranged in a discrete manner, and a 
transmission circuit and a reception circuit sharing the plurality of 
antennas for transmitting/receiving signals. The reception circuit includes 
a reception signal separating unit separating a signal from a specific 
terminal among the plurality of terminals based on signals from the 

20 plurality of antennas, when a reception signal is received, and a reception 
propagation path estimating unit estimating a reception response vector of a 
propagation path from the specific terminal based on signals from the 
plurality of antennas, when reception of the reception signal is received. 
The transmission circuit includes a transmission propagation path 

25 estimating unit estimating a transmission response vector of a propagation 
path when a transmission signal is transmitted based on a result of 
estimation by the reception propagation path estimating unit, and a 
transmission directivity control unit updating the antenna directivity when 
the transmission signal is transmitted, based on a result of estimation by the 

30 transmission propagation path estimating unit. The transmission 

propagation path estimating unit includes: an extrapolation processing unit 
calculating the transmission response vector of a down link slot to the 
specific terminal, by an extrapolation process based on a plurality of 
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reception response vectors of up link slots of the specific terminal estimated 
by the reception propagation path estimating- unit; a Doppler frequency 
estimating unit estimating a Doppler frequency of the propagation path, a 
memory holding a plurality of parameters used for the extrapolation process 
5 determined in advance in accordance with the Doppler frequencies of the 

propagation path; and a selecting unit selecting a parameter corresponding 
to the estimated Doppler frequency among the held plurality of parameters 
and applying the selected parameter to the extrapolation process by the 
extrapolation processing unit. The Doppler frequency estimating unit 

10 includes: a correlation operating unit calculating a vector correlation value 
based on reception response vectors preceding and succeeding in time 
estimated by the reception propagation path estimating unit; and an 
estimating unit estimating the Doppler frequency corresponding to the 
vector correlation value calculated by the correlation operating unit, based 

15 on the correspondence between vector correlation values and Doppler 
frequencies determined in advance experimentally. 

Therefore, not the reception signal itself but the correlation value of 
reception response vectors is found, and therefore, Doppler frequency of the 
propagation path for each separated specific terminal can be estimated 

20 correctly without the influence of the interfering component. Therefore, 
even when there is an estimation error in the reception response vector 
estimated for the up link, the transmission response vector of the down link 
can be estimated correctly, and hence, satisfactory control of transmission 
directivity can be realized. 

25 Preferably, the correlation operating unit includes a calculating unit 

calculating an instantaneous correlation value between reception response 
vectors preceding and succeeding in time and outputting the calculated 
value as the vector correlation value. 

Therefore, correlation between reception response vectors can be 

30 obtained at an arbitrary timing, and hence instantaneous Doppler frequency 
of the propagation path can be estimated correctly. Hence, more accurate 
estimation of transmission response vector becomes possible. 

More preferably, the correlation operating unit includes a 
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calculating unit calculating an instantaneous correlation value between the 
reception response vectors preceding and succeeding in time, and an 
averaging unit weight-averaging a past correlation value and a present 
correlation value calculated by the calculating unit with a prescribed weight 
5 coefficient, and outputting an obtained average value as said vector 
correlation value. 

Therefore, even when there is generated an error in the 
instantaneous Doppler frequency because of an abrupt fading, the Doppler 
frequency can be estimated correctly without the influence of such error, as 

10 the correlation values are averaged, and hence more accurate estimation of 
transmission response vector becomes possible. 

More preferably, the prescribed weight coefficient is set such that a 
weight for a past correlation value is large and a weight for a present 
correlation value is small. 

15 As the weight of the past correlation value is made large at the time 

of averaging, even when there is generated an error in the instantaneous 
Doppler frequency because of an abrupt fading, more accurate estimation of 
the Doppler frequency is possible, not influenced by the error, and hence 
more accurate estimation of transmission response vector becomes possible. 

20 More preferably, the correlation operating unit calculates a vector 

correlation value, based on a reception response vector of a present frame 
slot and a reception response vector of an immediately preceding frame slot. 

Therefore, different from the reference signal in the prior art, 
mutual correlation between reception response vectors can be obtained at 

25 any timing, and therefore, flexibility of operation to calculate the correlation 
value is increased. 

More preferably, the correlation operating unit calculates a vector 
correlation value based on a reception response vector of a present frame slot 
and a reception response vector of a most recent slot free of any reception 

30 error among past frame slots. 

Therefore, different from the reference signal in the prior art, 
mutual correlation between reception response vectors can be obtained at 
any timing, and therefore, flexibility of operation to calculate the correlation 
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value can be increased and, particularly, the influence of reception error can 
be eliminated. 

More preferably, the correlation operating unit calculates the vector 
correlation value based on a reception response vector of a former half and a 
5 reception response vector of a latter half of one same slot. 

Therefore, different from the reference signal in the prior art, 
mutual correlation of the reception response vectors can be obtained at any 
timing, and therefore, flexibility of operation to calculate the correlation 
value can be increased. 
10 The foregoing and other objects, features, aspects and advantages of 

the present invention will become more apparent from the following detailed 
description of the present invention when taken in conjunction with the 
accompanying drawings. 
BRIEF DESCRIPTION OF THE DRAWINGS 
15 Fig. 1 is a schematic block diagram showing a configuration of a 

radio equipment (radio base station) 1000 of a PDMA base station to which 
the present invention is applied. 

Fig. 2 is a flow chart illustrating an outline of the operation of the 
radio equipment (radio base station) 1000. 
20 Fig. 3 is a schematic diagram representing a basic operation of a 

transmission response vector estimator 32. 

Figs. 4A and 4B are illustrations representing a principle of 
estimating the transmission response vector of the present invention. 

Fig. 5 is an illustration representing the principle of determining 
25 extrapolation distance in accordance with an embodiment of the present 
invention. 

Fig. 6 is a flow chart representing an outline of the extrapolation 
process in accordance with an embodiment of the present invention. 

Fig. 7 is a schematic block diagram of a Doppler frequency 
30 estimating circuit in accordance with an embodiment of the present 
invention. 

Fig. 8 is a flow chart representing an example of a Doppler frequency 
estimating operation in accordance with an embodiment of the present 
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invention. 

Fig. 9 is a flow chart illustrating another example of the Doppler 
frequency estimating operation in accordance with an embodiment of the 
present invention. 

5 Fig. 10 is a flow chart representing a still further example of the 

Doppler frequency estimating operation in accordance with an embodiment 
of the present invention. 

Fig. 1 1 is a schematic block diagram of a Doppler frequency 
estimating circuit in accordance with another embodiment of the present 
10 invention. 

Fig. 12 is a flow chart representing an example of the Doppler 
frequency estimating operation in accordance with another embodiment of 
the present invention. 

Fig. 13 is a flow chart representing an outline of the extrapolation 
15 process in accordance with a still further embodiment of the present 
invention. 

Fig. 14 is a flow chart representing an outline of the extrapolation 
process in accordance with a still further embodiment of the present 
invention. 

20 Fig. 15 shows an arrangement of channels in various communication 

systems including frequency division multiple access, time division multiple 
access and PDMA. 

Fig. 16 is a schematic diagram representing a basic operation of an 
adaptive array radio base station. 
25 Fig. 17 is a schematic block diagram representing a configuration of 

an adaptive array radio base station. 

Fig. 18 is a schematic diagram representing a frame configuration of 
radio signals of a portable telephone. 

Fig. 19 is an illustration representing transmission/reception of 
30 radio signals between an adaptive array radio base station and a user. 
DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Fig. 1 is a schematic block diagram representing a radio equipment 
(radio base station) 1000 of a PDMA base station in accordance with an 
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embodiment of the present invention. 

In the configuration shown in Fig. 1, in order to distinguish a user 
PS1 from a user PS2, four antennas #1 to #4 are provided. The number of 
antennas, however, may be generally represented as N (N: natural number). 
5 In the transmission/reception system 1000 shown in Fig. 1, a 

reception unit SRI receiving signals from antennas #1 to #4 and separating 
a signal from a corresponding user (for example, user PSl), and a 
transmitting unit STl for transmitting a signal to user PSl are provided. 
Connection between antennas #1 to #4 and reception unit SRI or 

10 transmission unit STl is selectively switched by switches 10-1 to 10-4. 

More specifically, reception signals RXl(t), RX2(t), RX3(t) and 
RX4(t) received by respective antennas are input to reception unit SRI 
through corresponding switches 10-1, 10-2, 10-3 and 10-4, applied to a 
reception weight vector calculator 20 and a reception response vector 

15 calculator 22, as well as to one inputs of multipliers 12-1, 12-2, 12-3 and 12- 
4, respectively. 

To the other inputs of these multipliers, weight coefficients wrxll, 
wrx21, wrx31 and wrx41 for the reception signals of respective antennas are 
applied from reception weight vector calculator 20. The weight coefficients 
20 are calculated on real time basis by reception weight vector calculator 20, as 
in the prior art. 

Transmission unit STl includes: a transmission response vector 
estimator 32 receiving a reception response vector calculated by reception 
response vector calculator 22 and finding a transmission response vector by 

25 estimating a propagation path at the time of transmission, that is, 

estimating a virtual reception response vector at the time of transmission; a 
memory 34 communicating data with the transmission response vector 
estimator 32 and storing and holding data; a transmission weight vector 
calculator 30 calculating a transmission weight vector based on the result of 

30 estimation by transmission response vector estimator 32; and multipliers 
15-1, 15-2, 15-3 and 15-4 receiving at one inputs a transmission signal and 
at the other inputs, weight coefficients wtxll, wtx21, wtx31 and wtx41 from 
transmission weight vector calculator 30, respectively. Outputs from 
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multipliers 15-1, 15-2, 15-3, and 15-4 are applied through switches 10-1 to 
10-4 to antennas #1 to #4. 

Though not shown in Fig. 1, a configuration similar to that of 
reception unit SRI and transmission unit ST1 is provided for each user. 
5 [Principal of Operation of Adaptive Array] 

The operation of reception unit SRI will be briefly described in the 
following. 

The reception signals RXl(t), RX2(t), RX3(t) and RX4(t) received by 
the antennas are represented by the following equations. 

10 

RXi(t) = hiiSrxi(t)+hi 2 Srx 2 (t)+ni(t) ... (1) 

RX 2 (t) = h 2 iSrxi(t)+h 22 Srx 2 (t)+n 2 (t) ... (2) 

RXaCt) = h 3 iSrxi(t)+h32Srx 2 (t)+n 3 (t) ... (3) 

RX 4 (t) = h4iSrxi(t)+h42Srx2(t)+n 4 (t) ... (4) 

15 

Here, the signal RXj(t) represents a reception signal of the j-th (j = 1, 
2, 3, 4) antenna, while the signal Srxi(t) represents a signal transmitted 
from the i-th (i = 1, 2) user. 

Further, the coefficient hji represents a complex coefficient of the 
20 signal from the i-th user and received by the j-th antenna, while nj(t) 
represents noise included in the j-th reception signal. 

The equations (1) to (4) above in vector representation are as follows. 



X(t) = HiSr Xl (t)+H 2 Srx 2 (t)+N(t) ... (5) 

25 X(t) = [RXi(t) 5 RXz(t), RX 4 (t)] T -(6) 

Hi = [hii, h 2i , h4i] T , (i = 1, 2) ... (7) 

N(t) = [m(t), n 2 (t), n 4 (t)]T ... (8) 



In equations (6) to (8), [•••] T represents transposition of [ •••]. 
30 Here, X(t) represents an input signal vector, Hi represents a 

reception response vector of the i-th user, and N(t) represents a noise vector, 
respectively. 

The adaptive array antenna outputs a signal synthesized by 
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multiplying the input signals from respective antennas by weight 
coefficients wrxli to wrx4i, as a reception signal SRX(t) as shown in Fig. 1. 

Now, on the premise of the above described preparation, the 
operation of an adaptive array extracting the signal Srxl(t) transmitted 
from the first user is as follows. 

The output signal yl(t) of adaptive array 100 can be represented by 
the following equations, by a multiplication of input signal vector X(t) and 
the weight vector Wl. 



More specifically, the weight vector Wl is a vector that has, as 
components, weight coefficients wrxjl (j = 1, 2, 3, 4) to be multiplied by the 
j-th input signal RXj(t). 

Now, when we substitute input signal vector X(t) represented by the 
equation (5) for yl(t) represented by the equation (9), the following equation 
results. 

y 1 (t) = HiWi T Srxi(t) + H2W1T Srx 2 (t) + N(t) Wi T ... (1 1) 

When adaptive array 100 operates ideally here, weight vector Wl is 
sequentially controlled by weight vector control unit 1 1 to satisfy the 
following simultaneous equation, by a known method. 



When weight vector Wl is fully controlled to satisfy equations (12) 
and (13), the output signal yl(t) from adaptive array 100 is eventually given 
by the following equations. 



yi(t)=X(t)WiT 

Wi= [wrxii, wrx2i, wrx3i, wrx4i] T 




HiWi T = 1 
H 2 Wi T = 0 



... (12) 
... (13) 



yl(t) = Srxi(t) + Ni(t) 



- (14) 
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Ni(t) ni(t)wn + n 2 (t)w2i + ii3(t)w3i + n 4 (t)w4i 



- (15) 



Namely, as the output signal yl(t), the signal Srxl(t) transmitted by 
the first of the two users is obtained. 
5 [Summary of Operation of Radio Equipment 1000] 

Fig. 2 is a flow chart representing an outline of a basic operation of 
radio equipment 1000 as the base of the present invention. 

In radio equipment 1000, noting that the weight vector (weight 
coefficient vector) of the adaptive array can be represented uniquely by the 
10 reception response vector of each antenna element, the weight is indirectly 
estimated, by estimating time variation of the reception response vector. 

First, at reception unit SRI, propagation path of the reception signal 
is estimated, based on the reception signal (step S100). Estimation of the 
propagation path corresponds to finding an impulse response of the signal 
15 transmitted from the user, in accordance with equations (1) to (4). 

In other words, if the reception response vector HI can be estimated 
in equations (1) to (4), for example, the transmission path at the time of 
signal reception from user PS1 can be estimated. 

Thereafter, transmission response vector estimator 32 predicts 
20 propagation path at the time of transmission, that is, predicts reception 
response vector at the time of transmission, from the reception response 
vector at the time of reception (step S102). The predicted reception 
response vector corresponds to the transmission response vector at the time 
of transmission. 

25 Further, transmission weight vector calculator 30 calculates the 

transmission weight vector based on the predicted transmission response 
vector, and outputs the result to multipliers 15-1 to 15-4 (step S104). 
[Operation of Reception Response Vector Calculator 22] 
A basic operation as a base of the present invention of reception 
30 response vector calculator 22 shown in Fig. 1 will be described. 

First, assume that there are four antenna elements and two users 
communicate simultaneously. The signals output from reception circuit 
through respective antennas are represented by the equations (1) to (4) 
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above. 



At this time, the reception signals of the antennas represented by 
the equations (1) to (4) in vector representations are, again, given as 
equations (5) to (8). 



When the adaptive array is operating satisfactorily, signals from 
respective users are separated and extracted, and therefore, the signals 
Srxi(t) (i = 1, 2) all have the known values. 

At this time, utilizing the fact that the signal Srxi(t) is a known 
signal, it is possible to extract reception response vectors Hl = [hll, h21, 
h31, h41] and H2 = [hl2, h22, h32, h42] in the following manner. 

More specifically, by multiplying the reception signal by the known 
user signal, for example, the signal Srxl(t) from the first user, an ensemble 
average (time average) can be calculated in the following manner. 

E[X(t) • Srxi*(t)]=Hi-E[Srxi(t) • Srxi*(t)] 

+H2 • E[Srx 2 (t) • Srxi* (t) + E [N(t) • Srxi*(t)] ... (16) 

In equation (16), E[— ] represents the time average, and S*(t) 
represents a conjugate complex of S(t). When the time for averaging is 
sufficiently long, the average value can be given by the following. 

E[Srxl(t) • Srx*(t)] = 1 ... (17) 

E[Srx2(t) -Srxl*(t)] = 0 ... (18) 

E[N(t) • Srxl*(t)] = 0 ... (19) 

Here, the value of equation (18) becomes 0, as there is no correlation 
between the signal Srxl(t) and the signal Srx2(t). The value of the 



X(t) = HiSrxi(t) + H 2 Srx 2 (t) + N(t) 
X(t) = [RXi(t), RX2(t), RXn(t)] T 
Hi = [hu, ha, h ni ] T , (i = 1, 2) 
N(t) = [m(t), n 2 (t), n„(t)3 T 



... (5) 

- (6) 

- (7) 
... (8) 
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equation (19) becomes 0, as there is no correlation between the signal Srxl(t) 
and the noise signal N(t). 

Therefore, the ensemble average of equation (16) is, as a result, 
equal to the reception response vector HI, as shown below. 

E[X(t) • Srxi*(t)]=Hi ...(20) 

By the above described procedure, the reception response vector HI 
of the signal transmitted from the first user PSl can be estimated. 

In the similar manner, by an ensemble averaging operation of input 
signal vector X(t) and the signal Srx2(t), it is possible to estimate the 
reception response vector H2 of the signal transmitted from the second user 
PS2. 

The above described ensemble averaging is performed for a train of a 
prescribed number of data symbol at the head and a train of a prescribed 
number of data symbol at the end within one time slot at the time of 
reception. 

[Estimation of Transmission Response Vector] 

Fig. 3 is a schematic diagram representing a basic operation of the 
transmission response vector estimator 32 as a base of the present invention. 
Consider a 8-slot configuration in which 4 users are allocated for up and 
down links as PDMA burst. As to the slot configuration, the head 31 
symbols constitute the first training symbol train, succeeding 68 symbols 
constitute data symbol train and the last 3 1 symbols constitute the second 
training symbol train, for example. 

As described above, training symbol trains are provided at the head 
and at the end of the up link slot, and both reception response vectors are 
calculated using the algorithm of the reception response vector calculator 22 
described above. 

By extrapolation process (linear extrapolation), the reception 
response vector for the down link is estimated. 

More specifically, assuming that a value at an arbitrary time point t 
of an element of the reception response vector is f(t), based on the value f(t0) 
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at a time point tO of the head training symbol train of the up link slot and on 
the value f(tl) at a time point tl of the last training symbol train of the up 
link slot, the value f(t) at a time point t of the down link slot can be predicted 
in the following manner. 
5 f (t) = 

[f (tl) - f (tO)] / (tl - tO) x (t - tO) + f (tO) 

More specifically, the training symbol trains are provided at the 
head and at the end of the up link slot and linear extrapolation is performed 
in the above described example. A training symbol train may be 

10 additionally provided at a central portion of the up link slot, and the value 
f(t) at the time point t may be estimated by secondary extrapolation, from 
values of three p oints within the up link slot of the reception resp onse vector. 
Alternatively, higher order of extrapolation becomes possible by increasing 
the number of positions at which training symbol trains are provided within 

15 the up link slot. 

The present invention relates to an improvement of the method of 
estimating a reception (transmission) response vector for the down link by 
such an extrapolation process as well as to a Doppler frequency estimation 
circuit necessary for the improvement, details of which will be described 

20 later. First, determination of the transmission weight vector will be 
described. 

[Determination of Transmission Weight Vector] 
When the estimation value of the reception response vector at the 
time of transmission is calculated in the above described manner, the 
25 transmission weight vector can be found by any of the following three 
methods. 

i) Method Utilizing Orthogonalization 

Consider a weight vector W (1) (i) = [wtxll, wtxl2, wtxl3, wtxl4] at 
a time point t = iT (i: natural number, T: unit time interval) of a user PSl. 
30 In order to place nulls to user PS2, the following conditions should be 
satisfied. 

Assume that a predicted propagation path (reception response 
vector) for the user PS2 is V(2) (i) = [hi' (2) (i), h2' (2) (i), h3' (2) (i), h4' (2) (i)]. 
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Here, hp' (q) (i) represents a predicted value at a time point i of the reception 
response vector for the p-th antenna of the q-th user. It is assumed that the 
propagation path V(l) (i) for the user PSl has already been predicted in the 
similar manner. 

Here, W(l) (i) is determined to satisfy that W(l) (i) TV(2) (i) = 0. 
The following conditions cl) and c2) are imposed as binding conditions. 

cl) W(l) (i) TV(1) (i) = g (constant value) 

c2) || W(l) (i) || is a minimum value 

The condition c2) corresponds to minimization of the transmission 

power. 

ii) Method Utilizing Spurious Correlation Matrix 
Here, the adaptive array includes a number of antenna elements 
and a portion controlling the weight values of respective elements, as 
already described. Generally, when we represent the input vector of an 
antenna as X(t) and weight vector as W, optimal weight Wopt can be given 
by the following equation (Wiener solution) when the weight vector is 
controlled to minimize means square deviation between the output Y(t) = 
WTX(t) and the reference signal d(t) (MMSE standard: least square method 
standard). 



Here, YT represents transposition of Y, Y* represents complex area 
of Y, and E[Y] represents an ensemble average. By this weight value, the 
adaptive array comes to generate an array pattern to suppress unnecessary 
interfering wave. 

In the method utilizing the spurious correlation matrix, the equation 
(2 1) above is calculated in accordance with the spurious correlation matrix 




... (21) 



Here, the following relations must be satisfied. 



R X x = E[X*(t)xT(t)] 
rxd = E[x*(t)d(t)] 



... (22) 
... (23) 
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which will be described below. 

More specifically, the weight vector W(k) (i) for the user k is 
calculated, using the estimated complex reception signal coefficient h' (k)n 
(i). Assuming that the array response vector of the k-th user is V(k) (i), it 
5 can be given by the following equation, as already described. 

K ( * ) (0 = [*- 1 ( * ) (iX* 2 ( * ) W...*V ( * ) (o] ' -(24) 

At this time, the autocorrelation matrix Rxx (i) of the virtual 
10 reception signal at t = iT can be given by the following equation, using V(k) 
(i). 

^(0 = E^ ( " r (0^ )r (0+M ... (25) 

k=i 



15 Here, N is a virtual noise term, added as Rxx (i) should be an integer. 

In the calculation for the present invention, N is N = 1.0 x 10 -5 , as an 
example. 

The correlation vector rxd (i) between the reception signal and the 
reference signal is given by the following equation. 

20 

r x d(i)=V*>*(i) ...(26) 

Thus, the down link weight at time t = iT can b calculated by the 
equations (21), (25) and (26). 
25 The inverse matrix operation of the equation (25) can be optimally 

calculated for the user k in accordance with the lemma of inverse matrix. 
Especially when there are two users, the weight can be calculated in 
accordance with the following simple equations. 

30 W< l) (i)=(j> 12 +N)V m *<j)-p l2 V™\i) ...(27) 
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W i2) (/) -{Pn + N) V™* (0 - P 21 V w * (0 



■ • (28) 



The method of calculating a weight vector when autocorrelation 
matrix is given in this manner is described, for example, in T. Ohgane, Y. 
5 Ogawa, and K. Itoh, Proc. VTC '97, vol. 2, pp. 725-729, May 1997 or Tanaka, 
Ogahne, Ogawa, Itoh, Technical Report of the Institute of Electronics, 
Information and Communication Engineers, vol. RCS98-117, pp. 103-108, 
Oct. 1998. 

iii) Method Directing a Beam to User PS 1 
10 Focusing on the point that the beam is directed to the user PS 1, what 

is necessary is to satisfy the following equation. 
W(l)(i)=V(l)(i)* 

By determining the weight vector at the time of transmission by any 
of the above described methods for transmission, it is possible to suppress 
15 degradation of transmission directivity of the down link generated by the 
time difference between the up an down links even in the TDD/PDMA 
method, assuming a dynamic Raleigh propagation path including angular 
spread. 

Figs. 4A and 4B are illustrations representing the principal of 
20 estimating the transmission response vector in accordance with the present 

invention. The ideal state shown in Fig. 4A is, basically, a simplification of 

the concept shown in Fig. 3. 

More specifically, based on the reception response vector 1 and 

reception response vector 2 that are reception response vectors at two points 
25 within the same slot of the up link calculated in step S 102 of Fig. 2 by 

reception response vector calculator 22 shown in Fig. 1, linear extrapolation 

is performed to the original transmission timing of the corresponding slot of 

the down link, whereby a correct transmission response vector of the down 

link can be estimated. 
30 Here, the "ideal state" of Fig.4A stands on the premise that there is 

no estimation error in reception response vectors 1 and 2. 



-29- 



However, when there is an error in reception response vector 2' 
because of an estimation error resulting from a noise or sampling error in 
reception response vector 2, as shown in "example including response vector 
estimation error" of Fig. 4B, and linear extrapolation is performed in the 
5 similar manner (with the same extrapolation distance) as in the "ideal state", 
based on the reception response vectors 1 and 2\ the transmission response 
vector at the transmission timing would be further deviated, resulting in an 
erroneous estimation of transmission response vector. 

Therefore, when the process of determining the transmission weight 

10 (step S 104 of Fig. 2) is performed by the transmission weight vector 

calculator 30 shown in Fig. 1 based on the erroneous transmission response 
vector, the resulting transmission weight would also be erroneous, causing 
an error in down link directivity, that is, a transmission error. Particularly, 
as the distance between the radio base station and the terminal is long, a 

15 slight error in directivity causes a significant transmission error. 

Therefore, in the present invention, assuming that there is an 
estimation error in the reception response vector of the up link, an 
appropriate parameter, especially the extrapolation distance, for 
extrapolation process is adjusted in accordance with the propagation 

20 environment of the propagation path, so as to estimate a correct 

transmission response vector of the down link and to realize correct 
transmission directivity. 

Fig. 5 is an illustration representing the principal of determining the 
extrapolation distance in accordance with an embodiment of the present 

25 invention. 

The propagation environment of the propagation path is represented 
by the variation of reception coefficient of the propagation path, that is, the 
degree of fading, for example. The degree of fading is represented as a 
physical amount by a so-called Doppler frequency (FD). 

30 The Doppler frequency FD in the propagation environment is 

estimated generally in the following manner, in accordance with the 
principal of estimating Doppler frequency in accordance with the present 
invention. Namely, a correlation value between two reception response 
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vectors preceding and succeeding in time with respect to each other of the 
reception signals for each user separated by the adaptive array processing is 
calculated. When there is no fading, the two reception response vectors are 
identical, and therefore the correlation value will be 1. When there is 
5 significant fading, the difference between the reception response vectors 

increases, and the correlation value decreases. Such a relation between the 
correlation values of the reception response vectors and the Doppler 
frequencies FD is calculated in advance experimentally and held in a 
memory, whereby Doppler frequency FD at a certain time point can be 

10 estimated by calculating the correlation value of the reception response 

vectors. The Doppler frequency estimating circuit in accordance with the 
present invention will be described in greater detail later. 

First, the principal of operation in accordance with an embodiment 
of the present invention for determining the extrapolation distance, focusing 

15 on the Doppler frequency representing the degree of fading, will be 
discussed. 

As already described, when the reception response vector 2 is 
deviated as represented by the reception response vector 2' because of 
estimation error, the longer the extrapolation distance, the larger the 

20 extrapolation error, and therefore, the vector will be further away from the 
original transmission reception response vector. 

Generally, the smaller the fading, namely, the lower the Doppler 
frequency FD, the smaller the variation of reception coefficient on the 
propagation path. Therefore, in such a case, the extrapolation distance is 

25 made shorter, so as to prevent extrapolation larger than the actual amount 
of variation. More specifically, when the Doppler frequency FD is low, a 
short distance extrapolation from reception response vector 2' to the point a 
represented by X is performed as shown in Fig. 5(A), and the transmission 
response vector at point a is estimated and regarded as the correct 

30 transmission response vector of the point b represented by X. 

By contrast, the larger the fading, namely, the higher the Doppler 
frequency FD, the larger the variation of the reception coefficient of the 
propagation path. Therefore, in such a case, sufficient extrapolation is 
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ensured by making longer the extrapolation distance. More specifically, 
when the Doppler frequency FD is high, extrapolation of a relatively long 
distance from the reception response vector 2' to the point c represented by X 
is performed as shown in Fig. 5(B), and the transmission response vector at 
5 the point c is estimated and regarded as a correct transmission response 
vector of the point d represented by X. 

Such a process is mainly executed by the transmission response 
vector estimator 32 shown in Fig. 1. Fig. 6 is a flow chart representing the 
extrapolation process, focusing on the Doppler frequency FD. 

10 Referring to Fig. 6, in step Si, the propagation path is estimated by 

reception response vector calculator 22 of Fig. 1 and, specifically, reception 
response vectors 1 and 2' of the up link are estimated. 

Next, in step S2, the degree of fading, that is, Doppler frequency FD 
is estimated, by the method which will be described in detail later. 

15 Next, in step S3, by the transmission response vector estimator 32 

shown in Fig. 1, the optimal extrapolation parameter, that is, extrapolation 
distance, in accordance with the Doppler frequency FD is determined by the 
method described with reference to Fig. 5. For this purpose, it is assumed 
that optimal extrapolation distances determined by previous measurement 

20 in accordance with high/low values of Doppler frequency FD have been held 
in advance in memory 34 of Fig. 1. 

Thereafter, in step S4, extrapolation process is performed by using 
the extrapolation parameter (extrapolation distance) determined in step S3 
above, and propagation path of the down link, that is, transmission response 

25 vector, is estimated. 

Finally, in step S5, the transmission weight is estimated based on 
the transmission response vector of the down link determined in step S4 
above, by the transmission weight vector calculator 30 shown in Fig. 1. 
As described above, in the present embodiment, as the optimal 

30 extrapolation distance is selected dependent on whether the Doppler 

frequency FD is high or low, correct transmission response vector can be 
estimated, even when there is an estimation error in the reception response 
vector of the up link. 
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Fig. 7 is a schematic block diagram representing a configuration of a 
Doppler frequency estimating circuit in accordance with an embodiment of 
the present invention, used for the estimation of Doppler frequency FD (step 
S2 of Fig. 6) described above. The principal of operation of the Doppler 
5 frequency estimating circuit in accordance with the embodiment of the 
present invention will be described with reference to Fig. 7. 

Referring to Fig. 7, a response vector estimating circuit 101 
corresponds to the reception unit SRI of Fig. 1 and, particularly, to the 
reception response vector calculator 22 of Fig. 1. 
10 Here, four slots for the up link and four slots for the down link such 

as shown in Figs. 3, 4 and 5, namely, a total of 8 slots, will be referred to as 
one frame. Such frames are continued time sequentially, realizing 
alternating communication of up and down links. 

The response vector estimating circuit 101 applies the method of 
15 ensemble averaging described above to the reception signal, so as to 

estimate the reception response vector at a present frame slot, and applies 
the result to a correlation operation and Doppler frequency estimating 
circuit 103 and to a memory 102. 

The correlation operation and Doppler frequency estimating circuit 
20 103 calculates correlation value of the reception response vector of the 

present frame slot estimated by response vector estimating circuit 101 and 
the reception response vector of the corresponding slot of the preceding 
frame held in memory 102. 

The correlation value a of the reception response vectors of two 
25 frames preceding and succeeding in time is defined by the following 
equation. 

a|=hih 2 H |/|hi| |h 2 i 

Here, h.2 H represents a complex conjugate of each component of h2 
which is further subjected to transposition. 
30 Further, hi (i = 1, 2) represents the reception response vector (ha, hi2, 

hi3, hi4> having phase amplitude information of every antenna element of 
frame i as its elements. 

Though it is difficult to find an accurate correspondence between the 
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correlation values calculated in this manner and the Doppler frequencies, 
rough correspondence can be found experimentally. For example, if the 
correlation value is within the range of 1 to 0.95, it is estimated that the 
Doppler frequency FD is FD = 0Hz. If the correlation value is in the range 
5 of 0.95 to 0.80, it is estimated that FD = 10Hz, for example. 

The rough correspondence between the reception response vector 
correlation values and the Doppler frequencies FD obtained experimentally 
in the above described manner is stored in advance in correlation operation 
and Doppler frequency estimating circuit 103, and based on the correlation 
10 value of the vectors calculated in accordance with the equation above, a 
corresponding Doppler frequency estimated value is selected and output 
from circuit 103. 

Such a process as shown in Fig. 7 is generally executed in a software 
manner, using a digital signal processor (DSP), for example. Fig. 8 is a flow 
15 chart representing the process performed by the circuit configuration shown 
in Fig. 7. In the process shown in Fig. 8, a vector correlation value between 
a reception response vector of the present frame slot and the reception 
response vector of the corresponding slot of the immediately preceding frame 
is calculated. 

20 First, in step S31, the reception response vector of the present frame 

slot is estimated. 

Thereafter, in step S32, whether the reception response vector 
estimated in step S31 is the first estimated reception response vector or not 
is determined. If it is the first estimated reception response vector, it is 
25 stored in the memory (memory 102 of Fig. 7) in step S35. 

If it is not the first estimated reception response vector, a correlation 
value CORE between the reception response vector of the corresponding slot 
of an immediately preceding frame held in the memory and the reception 
response vector of the corresponding slot of the present frame estimated in 
30 step S31 is calculated in step S33. 

In step S34, based on the correspondence between the vector 
correlation values and the Doppler frequencies FD calculated 
experimentally in advance and held in the above described manner, the 
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Doppler frequency FD corresponding to the calculated correlation value 
CORK, is estimated and output. 

In step S35, the reception response vector of the present frame slot 
estimated in step S31 is stored in the memory (memory 102 of Fig. 7). 
5 By repeating the above described steps S31 to S35, it is possible to 

continuously obtain instantaneous vector correlation values between 
reception response vectors of corresponding slots of two frames that precede 
and succeed in time, that is, continuous two frames. 

Fig. 9 is a flow chart showing a modification of the process shown in 

10 Fig. 8. 

In Fig. 9, when there is no reception error detected in step S41, the 
following operation is basically the same as in the example shown in Fig. 8, 
and therefore, description thereof will not be repeated. When a reception 
error is detected in step S41, the flow cannot proceed to the next step S31 

15 until detection of the error is stopped. 

When detection of the reception error is stopped, the operation from 
step S31 to step S35 is executed. Here, it is noted that the process of step 
S33' is different from the example of Fig. 8. More specifically, the response 
vector estimation of the slot including the reception error is excluded in step 

20 S41. Therefore, in step S33', correlation value CORR between the reception 
response vector of the most recent slot free of any reception error among the 
past frame slots stored in the memory and the reception response vector of 
the present frame slot estimated in step S31 is estimated. The following 
operations are the same as in the example of Fig. 8. 

25 In the example shown in Fig. 9, the influence of the reception error 

can be eliminated, enabling more accurate estimation of the Doppler 
frequency. 

Fig. 10 is a flow chart showing another process performed by the 
circuit configuration shown in Fig. 7. In the process shown in Fig. 10, a 
30 vector correlation value between a reception response vector of a former half 
and a reception response vector of a latter half of one slot is calculated. 

First, in step S41, presence/absence of a reception error is 
determined (step S41 may be omitted). If detection error is not detected, a 
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reception response vector of a leading edge of the present slot is estimated in 
step S42 and thereafter, a reception response vector at the trailing edge of 
the present slot is estimated in step S43. 

Thereafter, in step S44, a correlation value CORR between the 
5 reception response vectors at the leading and trailing edges mentioned 
above is calculated. 

In step S45, based on the correspondence between the vector 
correlation values and the Doppler frequencies FD found experimentally 
and held as described above, the Doppler frequency FD corresponding to the 
10 calculated correlation value CORR is estimated and output. 

As described above, in the embodiments shown in Figs. 8 to 10, 
different from the reference signal as an object of correlation value 
calculation in the prior art, the correlation value of reception response 
vectors preceding and succeeding in time can be obtained at any timing, and 
15 therefore, flexibility of operation to calculate the correlation value can be 
improved. 

Fig. 1 1 is a schematic block diagram representing a configuration of 
a Doppler frequency estimating circuit in accordance with another 
embodiment of the present invention. The circuit configuration shown in 
20 Fig. 1 1 differs from the circuit configuration shown in Fig. 7 in the following 
points. 

In the example of Fig. 7, correlation operation and Doppler 
frequency estimating circuit 103 finds an instantaneous correlation value 
between two reception response vectors preceding and succeeding in time, 

25 while in the example of Fig. 11, correlation operation and the Doppler 

frequency estimating circuit 105 averages calculated correlation values, so 
as to enable more accurate estimation of the Doppler frequency. 

More specifically, referring to Fig. 11, correlation operation and 
Doppler frequency estimating circuit 105 operates a correlation value 

30 between the reception response vector of present frame slot estimated by 
response vector estimating circuit 101 and the reception response vector of 
the corresponding slot of the last frame held in memory 102, and calculates 
an average value between the just operated correlation value and a part 



-36- 



average correlation value read from the memory 104 which has been 
calculated and stored in advance. 

Specifically, the present embodiment is configured such that an 
average value weighted by prescribed weight coefficients is obtained. For 
5 example, in circuit 105, an average value between the value obtained by 
multiplying the past average correlation value stored in memory 104 by a 
first weight coefficient and a value obtained by multiplying the present 
correlation value just calculated by a second weight coefficient is calculated. 

Here, a large coefficient such as 0.97, for example, is set as the first 
10 coefficient, and a small coefficient such as 0.03 is set as the second 
coefficient. 

Based on the average correlation value obtained by the weighted 
mean, the circuit 105 estimates and outputs the corresponding Doppler 
frequency, and stores the thus obtained present average correlation value in 
15 memory 104. 

Fig. 12 is a flow chart representing a process performed by the 
circuit configuration shown in Fig. 11. In the process shown in Fig. 12, a 
vector correlation value between the reception response vector of the present 
frame slot and the reception response vector of an immediately preceding 
20 frame slot is calculated, and weighted mean of the past correlation value 
and the present correlation value is further calculated. 

First, in step S51, the reception response vector of the present frame 
slot is estimated in step S51. 

Thereafter, in step S52, whether or not the reception response vector 
25 estimated in step S5 1 is the first estimated reception response vector is 

determined. If it is the first estimated reception response vector, it is stored 
in the memory (memory 102 of Fig. 11) in step S59. 

If it is not the first estimated reception response vector, a correlation 
value CORR NOW between the reception response vector of the immediately 
30 preceding frame slot held in the memory and the reception response vector 
of the present frame slot estimated in step S51 is calculated in step S53. 

Thereafter, in step S54, whether the correlation value calculated in 
step S53 is the first calculated correlation value or not is determined. If it 
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is the first calculated correlation value, an averaging operation in which 
correlation value CORR NOW is assumed to be the averaged value CORR 
AVE is executed in step S55. 

If it is not the first calculated correlation value, a weight-averaging 
5 process is performed in step S56, in which the past average correlation value 
CORR OLD read from the memory and the present correlation value CORR 
NOW are multiplied by weight coefficients a and 1 - a, respectively. 

In step S57, based on the correspondence between the average 
correlation values and the Doppler frequencies FD calculated 
10 experimentally and held in advance as described above, a Doppler frequency 
FD corresponding to the calculated average correlation value CORR AVE is 
estimated and output. 

In step S58, the present average correlation value CORR AVE 
calculated in step S55 or S56 is stored as the past average correlation value 
15 CORR OLD in the memory (memory 104 of Fig. 11). 

Further, in step S529, the reception response vector of the present 
frame slot estimated in step S51 is stored in the memory (memory 102 of Fig. 
11). 

By repeating the above described steps S51 to S59, it is possible to 
20 estimate the Doppler frequency based on the average vector correlation 

value. As the average of the vector correlation values is taken, even when 
an error is generated in the instantaneous Doppler frequency because of an 
abrupt fading, for example, the Doppler frequency can be estimated 
correctly without the influence of the error. Particularly, by setting the 
25 weight coefficients such that the weight for the past correlation value is 

heavier, even when there is an error generated in the instantaneous Doppler 
frequency because of an abrupt fading, the Doppler frequency can be 
estimated more accurately, not influence by the error. 

In the example shown in Figs. 1 1 and 12, averaging process based on 
30 the correlation between the reception response vectors of corresponding slots 
of the continuous frames is performed. However, as in the example shown 
in Figs. 9 and 10, averaging process based on the correlation between the 
reception response vector of the present frame slot and the reception 
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response vector of the most recent slot free of any reception error among past 
frame slots may be performed. Alternatively, averaging process based on 
the correlation between reception response vectors of the former half and the 
later half of one slot may be performed. 
5 As described above, when the Doppler frequency estimating circuit 

in accordance with the present invention is used, an optimal extrapolation 
distance can be selected dependent on whether the Doppler frequency FD for 
each terminal estimated correctly is high or low, and hence, even when there 
is an estimation error of the reception response vector of up link, the 

10 transmission response vector can be estimated correctly. 

The propagation environment of the propagation path can also be 
represented by a weight estimation error of the signal obtained from the 
adaptive array output. Such an error is represented by Mean Square Error 
(MSE) between the signal value obtained from the adaptive array output 

15 and the expected desired signal value. The smaller the MSE, the more 

ideal estimation of the weight vector has been done in the up link, and hence 
the higher is the accuracy of the adaptive array output signal. By contrast, 
the larger the MSE, the less optimal the estimated weight vector in the up 
link, and the lower is the accuracy of the adaptive array output signal. The 

20 method of calculating the MSE is well known, and therefore, description 
thereof will not be given here. 

Therefore, dependent on the magnitude of MSE, there will be an 
estimation error in the reception response vector of the up link, resulting in 
extrapolation error. 

25 Focusing on the MSE, a principal of operation of another 

embodiment in accordance with the present invention will be described in 
the following. 

Returning to Fig. 5(A), when MSE is large and there is a large 
estimation error in the reception response vector, the extrapolation distance 
30 is made shorter, because the longer extrapolation distance causes larger 

extrapolation error. More specifically, when MSE is large, a short distance 
extrapolation is performed from the reception response vector 2' to the point 
a represented by X, and the transmission response vector at the point a is 
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estimated and regarded as the correct transmission response vector of the 
point b represented by x. 

When MSE is small and there is a small estimation error in the 
reception response vector as shown in Fig. 5(B), the extrapolation distance is 
5 made relatively long, because extrapolation error is small even when 
extrapolation distance is long. More specifically, when MSE is small, 
extrapolation of a relatively long distance from reception response vector 2' 
to the point c represented by X is performed as shown in Fig. 5(B), the 
transmission response vector at the point c is estimated and regarded as the 
10 correct transmission response vector for the point d represented by X. 

Such a process is mainly executed by transmission response vector 
estimator 32 shown in Fig. 1. Fig. 13 is a flow chart representing the 
extrapolation process focused on the MSE. 

The flow chart shown in Fig. 13 is the same as the flow chart shown 
15 in Fig. 6 except the following point. 

More specifically, in step S12, MSE is calculated by the known 
method, and in step S13, by the transmission response vector estimator 32 
shown in Fig. 1, the optimal extrapolation parameter, that is, extrapolation 
distance in accordance with the MSE is determined by the method described 
20 with reference to Fig. 5. It is assumed that for this purpose, optimal 

extrapolation distances determined by previous measurement in accordance 
with the magnitude of MSE have been held in memory 34 of Fig. 1. 

The process other than this point is the same as described with 
reference to Fig. 6, and therefore, description thereof will not be repeated 
25 here. 

As described above, in the present embodiment, as the optimal 
extrapolation distance is selected dependent on the magnitude of MSE, 
transmission response vector can be estimated correctly even when there is 
an estimation error in the reception response vector of the up link. 
30 Fig. 14 is a flow chart showing the extrapolation process in 

accordance with a still further embodiment of the present invention. In the 
embodiment shown in Fig. 6, the optimal extrapolation distance is 
determined by the Doppler frequency FD, and in the embodiment shown in 
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Fig. 13, the optimal extrapolation distance is determined by the MSE. In 
the embodiment shown in Fig. 14, which will be described below, the optimal 
extrapolation distance is determined considering both the Doppler 
frequency FD and MSE. 
5 More specifically, in the present embodiment, basically, the 

extrapolation parameter (extrapolation distance) is temporarily determined 
based on the Doppler frequency FD and, thereafter, the extrapolation 
parameter is corrected based on the MSE to reach the final determination. 
Referring to Fig. 14, steps Si and S2 are the same as those of the 

10 flow chart shown in Fig. 6, and therefore, description thereof will not be 
repeated. When the Doppler frequency FD is estimated in step S2, an 
optimal extrapolation distance is selected and temporarily determined, 
based on the correspondence between the Doppler frequencies FD and the 
extrapolation distances held in advance in memory 34. 

15 Thereafter, when MSE is estimated in step S12, the extrapolation 

distance is corrected in accordance with the magnitude of the MSE in step 
S22. For example, when the MSE is large, it is necessary to correct the " 
extrapolation distance to be shorter, and hence correction is done with a 
coefficient of X < 1. When the MSE is small, it is necessary to correct the 

20 extrapolation distance to be longer, and hence correction is done with the 
coefficient of X > 1. It is assumed that such coefficients have been found 
experimentally and held in memory 34. 

When the extrapolation distance is finally determined in step S22 in 
this manner, the transmission weight is estimated in steps SS4 and S5 that 

25 have been already described with reference to Fig. 6. 

As described above, in the present invention, a temporary 
extrapolation distance is selected in accordance with the Doppler frequency 
FD that has much influence on the extrapolation error, and the 
extrapolation distance is corrected by the MSE. Therefore, more accurate 

30 estimation of the transmission response vector is possible. 

The correspondence between the Doppler frequency or the MSE 
representing the propagation environment and the extrapolation distance 
differ equipment by equipment because of individual difference of the radio 
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equipment box. Therefore, generally, the correspondence is measured in 
advance for every radio equipment and determined. If it is assumed that 
there is only a slight individual difference, a common correspondence may 
be utilized for a plurality of radio equipments for the overall system. 

The extrapolation parameter is not limited to the extrapolation 
distance described above, and it may be other parameter including 
inclination of extrapolation. 

As described above, according to the present invention, even when 
there is an estimation error in the reception response vector estimated for 
the up link, it is possible to correctly estimate the transmission response 
vector of the down link by selecting a parameter used for the extrapolation 
process in accordance with the propagation environment, and hence, 
satisfactory transmission directivity can be realized. 

Further, according to the present invention, the Doppler frequency is 
estimated based not on the reception signal itself but on the correlation 
value of reception response vectors. Therefore, more accurate estimation of 
the Doppler frequency is possible without the influence of interfering 
component of the reception signal, while flexibility of operating process is 
increased. 

Although the present invention has been described and illustrated in 
detail, it is clearly understood that the same is by way of illustration and 
example only and is not to be taken by way of limitation, the spirit and scope 
of the present invention being limited only by the terms of the appended 
claims. 
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